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INTRODUCTION 


REPORT  NO  NCRE/R630  ADDENDUM 


This  addendum  contains  a precis  of  R630,  plus  additional 
computing  times,  control-cards  and  permissible  problem  size, 
regarded  as  superseding  that  document  for  the  conversant  SPAN 


information  on 
It  may  be 
user. 


The  loose-leaf  binding  is  intended  to  allow  for  easy  updating  of  the 
specification  as  and  when  program  developments  are  implemented.  The  data- 
card specification  given  in  R630  is  no  longer  valid  and  must  not  be  used. 


NCRE  PROGRAM  NO  118 


UNLIMITED 


<*t> 


Title : 


Computer: 
Language : 
Purpose : 


Method : 


^ IS 

"SPAN"  /i  a computer  program  for  static  and  dynamic  analysis 
of  stiffened  plates  and  grillages. 

CDC  6600  (SIA,  London). 

FORTRAN. 

The  analysis  of  plane  grillages  and  panels  under  normal  and 
in-plane  static  and  dynamic  loads.  Any  of  the  following 
operations  may  be  carried  out : 

(i)  evaluation  of  deflections  and  stresses  under  distrib- 
uted or  concentrated  lateral  loads. 

(ii)  evaluation  of  deflections  and  stresses  under  lateral 
loads  combined  with  in-plane  compressive  and/or  shear 
loads  (including  representation  of  initial  distortion). 

(iii)  evaluation  of  natural  frequencies  and  modes  of  vibra- 
tion. 

(iv)  analysis  of  deflections  and  stresses  under  steady- 
state  (harmonic)  dynamic  excitation. 

(v)  analysis  of  dynamic  response  to  general  time-dependent 
forces  (including  impulsive  loads). 

(vi)  evaluation  of  elastic  buckling  loads  and  modes. 


In  each  case,  the  program  is  applicable  to: 

(i)  flat  panels  of  unstiffened  isotropic  or  anisotropic 
plating  represented  by  plate  elements  only. 

(ii)  flat  stiffened  panels  represented  by  a combination 
of  beam  and  plate  elements. 

(iii)  grillages  represented  by  beam  elements  only. 

Full  details  of  the  methods  employed  are  given  in  NCRE 
report  R630.  For  present  purposes  it  is  sufficient  to 
say  that  the  program  is  based  on  the  finite  element  method; 
the  elements  themselves  being  represented  by  two-node 
Bernculli-Euler  beams,  with  shear  and  torsion  stiffness 
added,  and  triangular  three-node  Zienkiewicz  plates. 
Geometrical  representation  of  the  structure  is  aligned  to 
the  standard  concept  of  discrete  structural  elements  con- 
nected together  at  node  points , to  which  constraint  and 
load  conditions  are  also  referred. 


General  Notes : A data  generation  facility  is  incorporated  in  the  program. 

This  may  be  used  either  to  generate  a full  data  set  in- 
core  and  then  proceed  directly  with  the  calculation,  or 
to  punch  a full  data  set  on  cards,  perhaps  for  submission 
to  the  computer  after  some  minor  alteration.  The  second 


L 
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UIIUMITCn 


Idealisation : 


UNLIMITED 

usage  is  particularly  helpful  when  only  slight  departures 
from  the  permitted  regular  structure  and  load  config- 
urations occur. 


Data  specification  is  given  below  in  three  sections, 
"General",  "General  using  Datagen"  and  "Specific  Data 
for  Each  Calculation  Type".  Although  cards  must  be 
prepared  in  that  order,  the  last  section  should  be  per- 
used first  to  ascertain  the  calculation  type  number  and 
to  give  background  on  the  total  data  requirement . 


The  structure  is  reduced  to  a set  of  node  points  with 
interconnecting  beam  and  plate  elements,  as  appropriate. 

The  nodes  are  numbered  consecutively  from  1 to  NN,  and 
a right-handed  cartesian  axis  system  should  be  assumed , 
with  rotations  and  applied  moments  taken  as  clockwise 
positive  when  looking  in  the  co-ordinate  directions.  It 
is  possible  to  change  to  a local  axis  system  by  specifying 
nodal  transformations  from  the  initial  system.  If  this 
option  is  used,  applied  forces  and  restraints,  and  cal- 
culated deflections  are  all  referred  to  the  local  system. 

Note  that  the  structure  to  be  analysed  is  assumed  to  lie 
in  the  x-y  plane  and  to  possess  three  degrees  of  freedom 
at  each  node  point.  These  are  deflection  in  the  z direction 
and  rotations  about  the  x and  the  y directions,  taken  in  that 
order. 


Out-of-plane  loads  are  applied  as  either  point  forces  and 
moments  acting  at  nodes,  line-loads  (uniformly  distributed 
forces)  acting  on  beams,  or  normal  pressure  acting  on 
plates . In-plane  loads  take  the  form  of  an  axial  force 
for  beams,  and  two  direct  stresses  and  a shear  stress  for 
plates.  Tension  is  assumed  positive  throughout. 


Before  attempting  to  idealise  the  structure,  the  section 
"Problem  Size"  should  be  consulted  for  information  on  the 
permissible  maximum  number  of  nodes,  number  of  elements 
etc. 


It  is  not  always  easy  to  allocate  suitable  stiffness  pro- 
perties to  the  individual  elements,  and  particular  mention 
is  made  here  of  the  situation  existing  when  a grillage  is 
comprised  of  a mixture  of  beams  and  plates.  Since  no 
eccentricity  is  allowed,  each  beam  inertia  must  be  calcu- 
lated assuming  a part  of  the  plate  breadth  to  act  with  the 
beam  as  an  effective  flange-.  The  plates  themselves  are 
specified  as  normal  (thus  accounting  for  local  plate- 
bending behaviour).  If  in-plane  stresses  are  present, the 
beam  axial  forces  are  calculated  on  the  basis  of  the  nominal 
stress  multiplied  by  the  enhanced  cross-sectional  area, 
derived  from  the  plate-beam  combination.  Again,  the  plates 
are  specified  as  normal.  At  present,  the  problem  of  deciding 
on  a value  to  assume  for  the  breadth  of  plate  acting  with  a 
given  beam  is  not  completely  resolved.  However,  the  follow- 
ing figures  are  suggested  as  an  approximate  guideline. 
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UNLIMITED 


ji 


UNLIMITED 


Nominal  Plate  Breadth  (=  b/t) 

Effective  Breadth 
Nominal  Breadth 

Plate  Thickness 

Compressive 

Conditions 

Tensile 

Conditions 

b/t  < 40 

1.0 

0.8 

b/t  = 40  -*■  50 

1.0  -*■  0.75 

0.8 

b/t  = 50  80 

0.75  -*•  0.5 

0.8 

b/t  > 80 

(40)  x (t) 

0.8 

Data  Cards : 


General 


Card  1 : 


Card  2 : 


Card  3 : 


Card  4: 


(Up  to  80  alphameric  characters) 

TITLE 

F0RMAT(I5) 

NFLAG  = 0 Unabridged  data  input  (E  format) 

= 1 Unabridged  data  input  (F  format) 

Datagen  facility  used.  See  later 

= -2) 

FORMAT ( 1015, E15. 8) if  NFLAG  = 0 
FORMAT ( 101 5, F. 10.0) if  NFLAG  = 1 
NC  Calculation  type  number  (1  -»  6) 

NN  Number  of  Nodes 

NB  Number  of  beam  elements 

NP  Number  of  plate  elements 

NL  Number  of  lateral  load  conditions 

NNR  Number  of  restrained  nodes 

NNL  Number  of  loaded  nodes 

NBSEC  Number  of  beam  sections 

NPSEC  Number  of  plate  sections 

N00T  Output  parameter  (=  99  if  diagnostic  printout 

required) 

(=  0 otherwise) 

VI  Axial  load  parameter  (constant  multiplier 
acting  on  beam  axial  forces,  normally  set 
equal  to  one) 

omit  if  NNR  = 0 

F0RMAT( 4E15 . 8 )if  NFLAG  = 0 

FORMAT (4F10.0) if  NFLAG  = 1 
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UNLIMITED 


Card  8 ; 


Card  9*: 


omit  cards  8,  9 and  10  if  NP  = 0 
F0RMAT( 2E15 . 8 ) 

C0-0RD(1,I)  X Co-ord  of  node 

C0-0RD(2,I)  Y Co-ord  of  node 

- NN  cards,  ie  one  for  each  node. 
Specify  sequentially,  commencing 
at  one . 

FORMAT (E10.2,3I5,5E10.2) 

TH(I)  - Plate  thickness 

(-ve  if  stiffness  identical  to  previous 
plate) 


NXYC0(1 ,1 ) ) 
NXYC0( 2,1)  | 
NXYC0( 3,1)  ) 
XYRES(1,I)  ) 
XYRES( 2,1)  ) 


TRNl(I) 
TRN2( I ) 
TRN3( I ) 


) 


Node  numbers  of  plate  element 
vertices,  listed  anti-clockwise  in 
the  positive  Z direction 

Position  at  ) If  zero,  centroidal 
which  stresses)  values  assumed  by 
are  required  ) SPAN 

Angles  in  radians  between  local  nodal 
and  global  XY  planes.  Measure  from 
local  to  global,  clockwise  positive 
about  Z,  RH  axes  assumed  (see  Figure  1) 


Card  10*:  omit  if  TH  < 0 

FORMAT ( 3E1 5 . 8/6E10 . 2 ) 
SIG(l)  membrane  stress  o_ 

SIG(2)  membrane  stress  a, 

SIG(3)  membrane  stress  o 


xx 


yy 


xy 


Tensile  positive 


RPEC(1  ->-6)  plate  elastic  rigidities 
eg  isotropic  case:  h = TH,  E = Young's  Modulus, 
v = Poisson's  Ratio 
RPEC(l)  RPEC( 2 ) RPEC ( 4 ) 

RPEC( 3 ) RPEC(5)= 

SYM  RPEC ( 6 ) 


Eh 

12(l-v2) 


vEh" 


12(l-v  ) 
3 


Eh" 


12(l-v  ) 


SYM 


*(N0TE:  repeat  cards  9 and  10  NP  times). 


Gh" 

12 


8 - 

UNLIMITED 


UNLIMITED 

Card  11**:  omit  if  NB  = 0 


FORMAT (4E1 5 . 8/3E15 . 8 )if  NFLAG  = 0 
FORMAT (7F10.0) if  NFLAG  = 1 

BEAM  ELEMENT  NUMBER  (-VE  if  stiffness  and  distrib- 
uted load  identical  to  previous 
beam,  otherwise  no  limitations) 

NODE  NUMBER  at  x = 0 (local  beam  axis) 


0 at  x = 0 (angle  in  degrees  between  local  beam  axis 
and  local  nodal  axis . Measure  from 
local  nodal  to  local  beam  axis  - clock- 
wise positive  about  z , RH  axes  assumed 
(see  Figure  1)) 


NODE  NUMBER  at  x 


b (this  must  be  the  higher  of 
the  2 beam  node  numbers) 


0 at  x = b 


< = 
( 

(_ 

r 

OUTPUT  PARAMETERS  £ 

( = 
( 

(_ 

C 


0,  No  beam  stresses  or  forces 
output 

2 , Beam  stresses , forces  and 
displacement  output  for 
both  ends 

3 , Ditto  for  both  ends  plus 
mid-position 

n.  Ditto  for  n equally  spaced 
positions 


AXIAL  FORCE  Tensile  positive 

Card  12**:  omit  if  NB  = 0 or  beam  el  No  specified  -VE 


FORMAT (4E1 5. 8) if .NFLAG  = 0 
FORMAT (4F10.0) if  NFLAG  = 1 


TYPE  NUMBER  (double  ended  = 1) 
SECTION  IDENTIFIER 


b (element  length) 

LATERAL  LINE  LOAD  PARAMETER  (=  0,  no  lateral  line 

load 


(=  1,  otherwise 

Card  13**:  omit  if  NB  = 0 or -beam  el  No  specified  -VE,  or 

lateral  load  parameter  specified  zero 

FORMAT (4E1 5. 8) if  NFLAG  = 0 

FORMAT (8F10.0/8F10.0) if  NFLAG  = 1 


w 


NL 


) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 


Line  load  intensities  for  each  of  NL 
load  conditions 


A 


UNLIMITED 


UNLIMITED 


Data  Cards : 


**( NOTES:  1.  Cards  11  -*■  13  repeated  NB  times 

2.  If  NFLAG  = 0,  card  11  is  split  into 
2 cards , 11a  and  lib  - see  format 

3.  If  NFLAG  = 1,  card  13  is  split  into 
13a  and  13b  - see  format). 

NOW  REFER  TO  THE  SECTION  "SPECIFIC  DATA  FOR  EACH  CALCULATION 
TYPE"  FOR  DETAILS  OF  ANY  ADDITIONAL  CARDS  THAT  MAY  BE 
REQUIRED. 

General  Using  Datagen 


Card  1 : 


Card  2 : 


Card  3 : 


Card  4 : 


Card  5: 


(Up  to  80  alphameric  characters) 

TITLE 

FORMAT  (15) 

NFLAG  = -1  Full  data  set  generated  within  the 
computer  and  solution  attempted 

= -2  Full  data  set  punched  on  cards . No 
solut ion 

FORMAT ( 6110 ,F10.0 ) 

ND  Data  type  parameter  (1  or  2,  see  Figure  2) 

NC  Calculation  type  number  (1  ->•  6) 

NL  Number  of  lateral  load  conditions 

NOOT  Output  parameter  (=  99  if  diagnostic  printout 

required) 

(=  0 otherwise) 


) = 2,  Beam  stresses, 

, forces  and  dis- 


NOUTX  Output  for  x beams  ) (=  0,  No  beam  stresses  or 
NOUTY  Output  for  y beams  ) £ forces  output 

,=  2,  Beam  stresses, 

, forces  and  dis- 

placement  output 
^ for  both  ends 

(=  3,  Ditto  for  both  ends 
( plus  mid-position 

n.  Ditto  for  n equally 
^ spaced  positions 

VI  Axial  load  parameter  (constant  multiplier 
acting  on  beam  axial  forces,  normally  set 
equal  to  one) 

FORMAT (3F10.0) 

RL  X direction  grillage  length)  complete  grillage, 
B Y direction  grillage  length)  before  symmetry 

H Thickness  of  plating  (uniform) 
omit  cards  5,  6 and  7 if  H = 0 
FORMAT (6F10.0) 

RPEC  (1  -+•  6)  plate  elastic  rigidities,  see  page  8 


UNLIMITED 


UNLIMITED 


Card  6:  FORMAT ( 3F10 .0) 

SIG(l)  o xx  ) 

) membrane  stress  (tensile 
SIG(2)  0 yy  ) positive) 

SIG(3)  oxy  ) 

Card  7 : omit  if  NL  = 0 

FORMAT (8F10.0) 

Q ( I ) uniform  lateral  pressure,  I = 1 -*•  NL 

Card  8:  F0RMAT(2I10) 

NBX  no  of  beams  parallel  to  x-axis 

NBY  no  of  beai..s  parallel  to  y^axis 

(NOTES:  1.  NBX  and  NBY  refer  to  the  complete 
grillage,  even  if  axes  of  symmetry 
are  to  be  assumed 

2.  Edge  beams  are  not  allowed  when 
using  the  data-generation  facility 

3.  If  NBX  and  NBY  are  both  -VE,  beams 
are  considered  to  be  imaginary  and 
only  plat  . data  is  generated). 

Card  9:  FORMAT ( 110, 3F10.0) 

I Side  1 edge  restraint  parameter  (see 
Figure  3) 

Z Elastic  spring  constant  ) 

0 Elastic  spring  constant  . __  , 

x re,  ) input  only  if 

. ) i - 4 

0y  Elastic  spring  constant  ' 

Card  10:  as  above,  but  for  side  2 

Card  11:  as  above,  but  for  side  3 

Card  12:  as  above,  but  for  side  4 

Card  13*:  omit  cards  13,  14  and  15  if  NL  = 0 

FORMAT (110) 

II  Concentrated  load  (=  0,  omit  cards  14  and  15 

data  flag  (=  1,  supply  data 

Card  14*:  FORMAT(IIO) 

t h 

K Number  of  loaded  nodes  for  i load  condition 
(=  0 if  no  loads  for  that  condition  - omit 
Card  15) 

Card  15*:  FORMAT  4(2I5,F10.0) 

NXYC0(1)  X beam  number  ) 

NXYC0(2)  Y beam  number  j 4 sets  per  card 

CNLF  Applied  nodal  force  ) 

(Specify  full  value  on  axis  of  symmetry) 

*( NOTES:  1.  Nodes  are  defined  by  beam  inter- 
sections, X beams  are  numbered  from 
0(Y  = 0 edge)  to  NBX+1.  Ditto  for 
Y beams. 


Unlimited 


Card  18: 


Card  19 : 


Card  20: 


UNLIMITED 

2.  forces  are  only  specified  for  nodes 
remaining  after  imposition  of 
symmetry,  if  any. 

3.  Cards  14  and  15  are  read  together, 
card  15  being  repeated  K times. 
Taking  cards  14  and  15  as  a group, 
there  are  NL  groups  of  data. 

TYPE  1 DATA 

X beam  data  (omit  if  NBX  =0) 


Card  16:  F0RMAT( 8F10 .0) 


DX(  I ) 

V ( M3 ) * 

V(M3+1)* 

V(M3+2)* 

V(M3+3) 

V(M3+4) 


Y direction  spacing  of  X beams 
Inertia  constant 
Torsion  constant 
Shear  area 

Distance  from  NA  to  outer  fibre 
Young's  modulus 


V(M3+5)  Shear  modulus 

RPX*  Axial  force  (tensile  positive;  if  zero 

punch  0) 

’•(NOTE:  Specify  full  values  on  axis  of  symmetry). 
Card  17:  omit  if  NL  = 0 or  NBX  and  NBY  both  -VE 

FORMAT (110) 

J = 0,  no  line  load 

= 1,  supply  NL  line  loads 


omit  if  NL  = 0 or  NBX  and  NBY  both  -VE  or  J = 0 
FORMAT (8F10.0) 

WX(JJ)  line  load  intensities,  supply  NL  values 
(specify  full  value  on  axis  of  symmetry) 

Y beam  data  (omit  if  NBY  =0) 

repeat  card  16  using  data  for  Y beams . If  Y beams 
have  same  section  as  X beams,  set  DY  -VE  (any  number) 
and  leave  rest  of  card  blank,  unless  last  X beam 
lies  on  an  axis  of  symmetry  (in  this  case  supply 
full  set  of  Y beam  data). 


Card  21-  ) rePeat  cards  17  aad  18  as  necessary 

TYPE  2 DATA 

X beam  data  (omit  if  NBX  =0) 

Card  16**;  FORMAT (8F10.0) 


DX(L)  Y direction  spacing  from  previous  X 

beam  to  present  X beam  (or  from  Y = 0 
for  first  beam) 

RPX(L)*  Axial  force  (tensile  positive;  if  zero 
punch  0) 


12  - 

UNLIMITED 


Data  Cards : 


UNLIMITED 


V(M3)*  Inertia  constant  (-1  if  similar  beam 
already  specified) 

V(M3+1)*  Torsion  constant  (or  section  number  of 
similar  beam  already  specified) 

(NOTES:  1.  Beam  section  numbers  are  allocated 
sequentially  from  1 for  all  beams 
with  positive  inertia  constant  . 

2.  Omit  the  following  card  16  entries 
if  the  inertia  constant  is  -1). 


V(M3+2)*  Shear  area 

V(M3+3)  Distance  from  NA  to  outer  fibre 

V(M3+4)  Young's  modulus 

V(M3+5)  Shear  modulus 

*(NOTE:  Specify  full  values  on  axis  of  symmetry). 
Card  17**:  omit  if  NL  = 0 or  NBX  and  NBY  both  -VE 


FORMAT ( 110) 

I = 0,  no  line  load 

= 1,  supply  NL  line  loads 

Card  18**:  omit  if  NL  = 0 or  NBX  and  NBY  both  -VE  or  I = 0 


FORMAT (8F10.0) 

WX(II)  line  load  intensities,  supply  NL  values 

(specify  full  values  on  axis  of  symmetry) 

**(N0TE:  Cards  16,  17  and  18  are  repeated  NSYMX 
times,  where  NSYMX  is  the  number  of 
beams  remaining  after  imposing  symmetry. 
Beams  lying  on  an  axis  of  symmetry  are 
included  in  this  figure). 

Y beam  data  (omit  if  NBY  =0) 

Card  19:  ) 

Card  20:  | repeat  cards  16,  17  and  18  as  necessary 

Card  21 : ) 


Specific  Data  for  Each  Calculation  Type 

Notation:  [k]  = matrix  of  bending  stiffness  coefficients 


w = 

[6]  = 
M = 


geometric  matrix  - contains  additional 
bending  stiffness  coefficients  arising 
from  in-planc  stresses 

deflection  vector 

static  load  vector 


X 

[M] 


eigenvalue 

diagonal  matrix  of  masses  and  mass  moments 
of  inertia 


[c] 

u> 

e 


diagonal  matrix  of  damping  coefficients 
vector  of  time  dependent  force  amplitudes 
angular  frequency  ) 
phase  angle  ) relatin8  to 


k 


13  - 


UNLIMITED 


Card  1: 


Card  2 : 


Card  3 : 


Card  4: 


Lateral  Load  Only 
Type  Number,  NC  = 1 
Governing  Eq  [k]  [fi]  = [r] 

No  additional  data  required 
Combined  Lateral  and  Axial  Load 
Type  Number,  NC  = 2 


Governing  Eq 


[k]  + [KJ  [6]  = [R] 


No  additional  data  required 
Natural  Frequencies  and  Mode  Shapes 
Type  Number,  NC  = 3 

Governing  Eq  [M  + [kg]  -A  Cm]]  [6]  = o 
FORMAT (15) 

IREAD  = 0 stiffness  matrix  read  from  cards 
(upper  triangle  only,  stored  in 
columns ) 

= 1 matrix  generated  by  SPAN 
omit  if  IREAD  = 1 


FORMAT( 15 ) 

NMAT  = order  of  matrix 
omit  if  IREAD  = 1 
FORMAT (8F10.0) 

K(I),  I : 1+M  stiffness  matrix,  where  M = (NMAT)  x 
(NMAT+l)/2,  (upper  half  only,  in  cols.  ) 

FORMAT( 215 ) 

NMODE  = number  of  vibration  modes  required  in 
result 


NDOF  = number  of  significant  degrees  of  freedom 

(relating  to  non-zero  masses) 

Card  5*:  FORMAT ( 15 ,F10.0) 

NDF(J)  = degree  of  freedom  number  (degrees  of 
freedom  re-numbered  after  imposing 
constraints  but  including  zero  mass 
(ie  non-significant)  values) 

SM(J)  = mass 

*(NOTE:  Repeat  Card  5 NDOF  times). 

Steady  State  Response 

Type  Number,  NC  = 4 

Governing  Eq  [m]  [fi]  + [c]  [«]  + jjx]  + QcQ]j  [«] 

= [F  cos(ut  + e)] 


UNLIMITED 


Card  1 : 


Card  2*: 


Card  1: 


Card  2: 


UNLIMITED 

FORMAT( 215 ,F10. 0) 


NDOF  = 

NDAMP  ( = 
( = 

w = 


number  of  significant  degrees  of  freedom 
(relating  to  non-zero  masses) 

0,  no  damping 

1,  damping 

forcing  function  frequency  (radians/sec) 


F0RMAT( 15 ,4F10 .0) 

NDF(I)  = significant  degree  of  freedom  number' 
(degrees  of  freedom  re-numbered  after 
imposing  constraints  but  including 
zero  mass  (ie  non-significant)  values) 


SM( I ) = mass 

C( I)  = damping  coefficient 

F(I)  = exciting  force  amplitude 

PAF(I)  = phase  angle  relating  to  F 

*(N0TE:  Repeat  Card  2 NDOF  times). 

Transient  Response 


Type  Number,  NC  = 5 

Governing  Eq  [m]  [s]  + [c]  [<s]  t ^[k]  + 


FORMAT (4110) 


K2 

- 

number 

of  exciting  forces 

K3 

= 

number 

of  imposed  displacements 

K4 

= 

0 

K5 

( = 

0 

if  X 

is  zero  at  all  times 

( = 

1 

if  X 

is  non-zero  but  constant 

( = 

2 

if.  X 

is  time-dependent 

( = 

3 

if  X 

is  constant  within  each 

integration  range 

(NOTE:  X is  a .scale  factor  acting  on  K , the 
geometric  stiffness) 


FORMAT (8F10.0)  as  many  cards  as  necessary 

List  NMAT  masses  followed  by  NMAT  damping 
coefficients,  where  NMAT  is  the  order  of  the 
stiffness  matrix  after  constraining.  Do  not 
start  a new  card  for  the  damping  coefficients 
unless  the  previous  card  is  full. 


(NOTE:  If  zero  masses  are  given,  the  corresponding 
degrees  of  freedom  are  condensed  out  of  the 
solution  and  re-inserted  later  on,  in  the 
results  stage). 


■■1 


15 


UNLIMITED 


Card  3: 


Card  4: 


Card  5 : 

Card  6 : 


Card  7 : 


Card  8 : 


omit  if  K^  = 0 

FORMAT (F10.0) 

X 

omit  if  Kg  = 0 
F0RMAT(8F10.0) 

Degree  of  freedom  numbers  for  imposed  displacements 
FORMAT (110) 

NIR  = number  of  integration  ranges 
FORMAT ( 8F10 . 0 ) 

Imposed  displacements  vector  followed  by  Initial 
velocity  vector,  of  length  NMAT  in  each  case 

(do  not  start  velocity  vector  on  a new  card  unless 
previous  one  full.  If  vectors  null , punch  zero's) 

omit  if  Kg  = 0 
F0RMAT(4(I10,F10.0) ) 

Degree  of  freedom  number  for  forces , and  force 
magnitudes 

omit  if  Kg  = 0 
F0RMAT(4( 110  ,F10.0) ) 


Degree  of  freedom  number  for  imposed  displacements , 
and  displacement  magnitudes 

Card  9:  FORMAT ( 110, F10. 0,3110) 

Number  of  time  seeps 

time  step 

Number  of  time  steps  between  each  printout 


6 control  parameter : 

(=  0,  compute  new  F 

(=  1,  compute  new  reduced  F K 

(=  2,  compute  new  reduced  F K M C 

K (=  0,  if  exciting  forces  and  imposed 

displacements  constant 

(=  1,  if  exciting  forces  and  imposed 

displacements  time-dependent 

Card  10*:  FORMAT (F10.0) 


X , needed  if  K_  = 2 and  K.  i 0 
o b 

Card  11*:  F0RMAT(4(I10,F10.0) ) 

Degree  of  freedom  number  for  forces , and  force 
amplitudes;  needed  only  if  Kg  and  Kg  non-zero. 
Specify  forces  Kg  times 


UNLIMITED 


Card  12*: 


Card  13: 


Card  1: 


Card  2: 


Card  3: 


Card  4: 


UNLIMITED 

F0RMAT(4(I10,F10.0) ) 

Degree  of  freedom  number  for  displacements , and 
displacement  amplitudes ; needed  only  if  and 
non-zero,  times 

*(N0TE:  Cards  10,  11  and  12  are  repeated  NTS 
times) • 

FORMAT (F10.0) 

X,  required  for  second  and  subsequent  integration 
ranges  (needed  only  if  K = 3) 

Repeat  data  from  card  7 for  next  integration  range , 
as  necessary 

Linear  Buckling** 


Type  Number,  NC  = 6 
Governing  Eq  [W  * 
F0RMAT( 315 ,2F10.0) 


eg]  w -- 


NMODE  Number  of  buckling  modes  required  in 
results.  If  left  unspecified,  value 
set  by  program  = 4 

IREAD  Option  to  read  in  matrices  K and  K . 

Specify  zero  if  not  required,  1 otherwise 

ICON  Number  of  degrees  of  freedom  to  be  con- 
densed out  of  the  solution 

GP  ) Lower  and  upper  limits  for  the  required 
V buckling  loads*  (order  immaterial). 

^ If  left  unspecified,  GP  = 0.001,  GQ  = 1000 

*(N0TE:  Use  the  reciprocal  of  the  anticipated 
loads  for  the  purpose  of  specifying 
limits) . 

omit  if  ICON  = 0 

FORMAT (1615) 

N(I)I=l->-  ICON.  Degree  of  freedom  numbers 

(refers  to  re-numbered  set, 
after  imposing  constraints). 
Repeat  card  as  necessary. 

omit  if  IREAD  =’0 
F0RMAT(I5) 

NMAT  Order  of  matrices  to  be  read  in 
omit  if  IREAD  = 0 
FORMAT (2F10.0) 

K(  I ) , K„(I)  1 = 1-*-  NMAT  (NMAT  + l)/2 

O 

Conventional  and  geometric  stiffness  matrices. 
Supply  upper  half  only  (in  columns) 

**(N0TE:  It  should  be  remembered  that  the  buckling 
load,  calculated  and  printed  by  the 
program,  is  a uniform  multiplier  acting 
on  the  beam  axial  forces  and  plate 
stresses.  The  actual  buckling  load  is 
arrived  at  by  multiplying  the  printed 
buckling  load  by  the  values  allotted  to 
the  input  parameters  "SIG(1,2,3)", 

"AXIAL  FORCE",  or  "RPX",  as  appropriate). 

- 17  - 

UNLIMITED 


UNLIMITED 


Output : 


General 


For  all  calculation  types,  the  title  is  printed  out  followed 
by  the  value  of  the  control  parameter  "NFLAG".  If  NFLAG  = -1 
(ie  the  datagen  facility  i"  used)  a list  of  generated  input 
data  is  printed,  followed  by  an  array  containing  the  degree 
of  freedom  numbers  eliminated  by  constraint  conditions  (the 
NDE  array).  If  NFLAG  = -2  a list  of  generated  input  data  is 
printed,  and  an  unabridged  card  deck  is  punched  in  E format. 
The  run  is  then  terminated. 

Diagnostic  output  may  be  obtained  by  setting  the  parameter 
"NOOT"  equal  to  99. 

Lateral  Load  Only,  and  Combined  Lateral  and  Axial  Load 

After  the  general  heading  "NODAL  DISPLACEMENTS  AND  ROTATIONS" 
the  following  results  are  printed  for  each  lateral  load 
condition: 

a.  Deflections  and  rotations  for  each  node. 

b.  Beam  element  displacements,  moments,  forces  and  bending 
stresses  at  each  stress  point  (specified  in  the  input 
data) . 

c.  Plate  element  moments,  computed  at  centroids  when  the 
datagen  facility  is  used. 

Natural  Frequencies  and  Mode  Shapes 

Frequencies  (radians/sec)  and  mode  shapes  are  printed  out  for 
the  lowest  NMODE  resonant  frequencies. 

Steady  State  Response 

The  following  results  are  printed: 

a.  Displacement  amplitudes  and  phase  angles  (radians)  for 
each  degree  of  freedom  (remaining  after  imposing  con- 
straints). 

b.  Beam  element  displacement,  moment,  force  and  bending 
stress  amplitudes,  together  with  corresponding  phase 
angles,  at  each  stress  point. 

c.  Plate  element  moment  amplitudes  and  phase  angles,  com- 
puted at  centroids  when  the  datagen  facility  is  used. 

Transient  Response 

The  following  results  are  printed : 

a.  List  of  options  selected. 

b.  Tables  giving  bending  stiffnesses  (K  ),  geometric 
stiffnesses  (K  ) and  modified  stiffness  matrices. 

c.  For  each  time  step,  the  displacement,  velocity  and 
acceleration  for  each  degree  of  freedom. 


i8UNUM!I£fi 


UNLIMITED 


d.  If  the  calculation  includes  more  than  one  integration 
range,  'b'  and  'c'  are  repeated. 

Linear  Buckling 

Buckling  loads  and  mode  shapes  are  printed  out  for  the  lowest 
NMODE  eigenvalues.  The  modal  vector  contains  deflections, 
listed  in  rows , for  those  degrees  of  freedom  that  remain 
after  imposing  constraints  and  carrying  out  condensations 
(if  any). 

Problem  Size:  General 


When  using  the  datagen  facility  the  following  limitations 
apply: 

1.  Maximum  no  of  nodes  = 91. 

2.  Maximum  no  of  loadcases  = 10. 

3.  Total  no  of  X and  Y beams,  taken  collectively,  should 
not  exceed  24. 

For  any  calculation: 

4.  Maximum  no  of  degrees  of  freedom  that  may  be  condensed 
out  of  the  solution  = 80. 

Calculation  Types  1 and  2 


1. 

Maximum 

no 

of 

nodes  = 91. 

2. 

Maximum 

no 

of 

beam  elements  = 

162. 

3. 

Maximum 

no 

of 

plate  elements  = 

150. 

4. 

Maximum 

no 

of 

restrained  nodes 

= 50. 

5. 

Maximum 

no 

of 

loadcases  (NL)  = 

2. 

For  smaller  calculations: 

6.  NL  = (5859  -(31xNB))/2,  8100/(27xNP) , or  768/(3xNN), 
whichever  is  least . 

Calculation  Types  3 and  6 

1.  Maximum  no  of  nodes  = 64. 

2.  Maximum  no  of  degrees  of  freedom,  remaining  after 
imposing  constraints  = 150. 

3.  Maximum  no  of  eigenvalues  (ie  buckling  loads  or 
resonant  frequencies)  = 5. 

4.  Maximum  no  of  beam  elements  = 162. 

5.  Maximum  no  of  plate  elements  = 150. 


- UNLIMITED 


UNLIMITED 


6.  (Type  3 
( 

(Type  6 
( 


Minimum  no  of  significant  degrees  of 
freedom  = 70. 

Maximum  no  of  insignificant  degrees 
of  freedom  (ICON)  = 80. 


For  smaller  calculations: 


7.  Maximum  no  of  eigenvalues  = 768  * (no  of  degrees  of 
freedom  remaining  after  imposing  constraints). 

8.  Type  3 - Minimum  no  of  significant  degrees  of  freedom 

(no  of  degrees  of  freedom  remaining  after 
imposing  constraints)  -80. 

Calculation  Type  4 

For  the  undamped  case,  type  3 limitations  apply. 


If  damping  is  specified: 


1.  Maximum  no  of  significant  degrees  of  freedom  = 75. 
The  minimum  value  is  zero. 


2.  In  other  respects,  type  3 limitations  apply. 


Calculation  Type  5 

1.  Maximum  no  of  degrees  of  freedom,  remaining  after 
imposing  constraints  = 36. 

2.  In  other  respects,  type  3 limitations  apply. 

Error  Messages : Various  data  checks  are  carried  out  by  SPAN,  and  if  the  data 
is  found  to  be  inadmissible  an  error  message  will  result. 
Advice  on  how  to  interpret  such  a message  is  given  in  NCRE 
Report  R630.  Generally,  experience  shows  that  FORMAT  errors 
or  missing  data  cards  are  by  far  the  most  probable  cause  of 
run  failure. 


Control : The  run  stream  is  constituted  as  follows: 

aaa ,CM305000 ,Taa . 

ACCOUNT ,aaaaaa ,aaaa ,aaa . 

VSN( OWNER  = 124900  ,T  = 2440) 

LABEL(T.R) 

C0PYBF(T,BSPAN) 

REWIND(BSPAN) 

LOAD(BSPAN) 

EXECUTE (SPAN) 

7/8/9 

data  cards 
6/7/8/9 


UNLIMITED 

Run  Times:  There  is  a considerable  variation  in  run  time,  for  a given 

problem  size,  depending  on  the  type  of  calculation  being 
executed.  The  following  table  gives  some  typical  figures: 


Calculation  Type  No  = 

D 

a 

3 

a 

5 

6 

No  of  Nodes  = 

9 

45 

8 

5 

4 

32 

No  of  Beams  = 

0 

52 

n 

4 

3 

0 

No  of  Plates  = 

8 

0 

0 

0 

42 

No  of  Degrees  of  Freedom 
in  Final  Solution  = 

12 

83 

9 

3 

63 

Total  No  of  Time  Steps  = 

/ 

/ 

/ 

750 

/ 

CPU  Time  (Secs)  = 

0.4 

3.0 

0.2 

26.9 

6.5 

1 


Programmed  and  Implemented  by  G C Mitchell  and  others 


NCRE/R630  ADD 
UNLIMITED  A 


UNLIMITED 


GRILLAGE  TYPES  FOR  DATAGEN  FACILITY 


TYPE  I -REGULAR  GR ILLAGE  (ORTHOGONAL J 


Y BEAMS 


X BEAMS 


NOTES-  I.  ALL  X BEAMS  HAVE  SAME  SECTION 
ALL  Y BEAMS 

2.  SPACINGS  DX  ARE  CONSTANT 
SPACINGS  DY  ii  " 


TYPE  2 — IRREGULAR  GRILLAGE  (ORTHOGONAL) 

NOTES  - I.  X AND  Y BEAMS  MUST  HAVE  CONSTANT  SECTIONAL 
PROPERTIES  ALONG  THEIR  LENGTHS 

2.  EACH  X AND  Y BEAM  MAY  HAVE  A DIFFERENT 
SECTION  FROM  ITS  NEIGHBOUR 

3.  SPACINGS  DX  AND  DY  MAY  VARY  FROM  BEAM 
TO  BEAM 


FIGURE  2 


UNLIMITED 


NCRE/R63Q  ADD. 


UNLIMITED 


GRILLAGE  EDGE  RESTRAINT  CONDITIONS 


CARDS  9,10,11,12: 

l*  |.  UNRESTRAINED 

2.  SIMPLY  SUPPORTED 

3.  EDGE  CLAMPED 

4.  EDGE  ELASTICALLY  RESTRAINED 

5.  PLANE  OF  SYMMETRY  AT  MID  SPAN 

6.  PLANE  OF  ANTI  - SYMMETRY  AT  MID  SPAN 

TYPES  5 AND  6 EDGE  RESTRAINTS  CAN  BE 
APPLIED  TO  SIDES  3 AND  4 ONLY 


SIDE  3 


FIGURE  3 


NCRE  / R63Q  ADD. 

UNLIMITED 


